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Introduction 



The theory of pseudodifferential operators (PDOs) is a powerful technique, which has 
many applications in analysis and mathematical physics. In the framework of this theory, 
one can effectively construct the inverse of an elliptic differential operator L on a closed 
manifold, its non-integer powers and even some more general functions of L. For operators 
with constant coefficients in IR n , this can be easily done by applying the Fourier transform. 
In a sense, the theory of PDOs extends the Fourier transform method to operators with 
variable coefficients and operators on manifolds at the expense of losing infinitely smooth 
contributions. This is normally acceptable for theoretical purposes and is useful for numerical 
analysis, since numerical methods for the determination of the smooth part are usually more 
stable. 

Traditionally, PDOs on manifolds are defined with the use of local coordinates. This leads 
to certain restrictions on operators under consideration, as all the definitions and results 
must be invariant with respect to transformations of coordinates. The main aim of this 
paper is to introduce the reader to a little known approach to the theory of PDOs that 
allows one to avoid this problem. 

The paper is constructed as follows. In Section [1] we recall some basic definitions and 
results of the classical theory of PDOs. Their detailed proofs (as well as other relevant 
statements and definitions) can be found, for instance, in [H2[ IShu[ ITa | ITr] . Section [2] gives 
a brief overview of some elementary concepts of differential geometry (see |KNj or any other 
textbook for details). In Sections [3] and H] we explain how to define PDOs without using 
local coordinates and quote some results from the paper [Sal] and the conference article 
[Sa2j. Section [5] contains new results on approximate spectral projections of the Laplacian 
obtained in the PhD thesis |McK] . Finally, in Section [6] we give a review of other related 
results and discuss possible developments in the field. 

Throughout the paper denotes the space of infinitely differentiable functions with 
compact supports, and T>' is the dual space of Schwartz distributions. Recall that, by the 
Schwartz theorem, for each operator A : Cq° i— > T>' there exists a distribution A(x, y) G T>' 
such that (Au, v) = (A(x,y) , u(y)v(x)) for all u,v G C£°. The distribution A(x, y) is called 
the Schwartz kernel of A. 
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I. PDOS: LOCAL DEFINITION AND BASIC PROPERTIES 

Let a(x, y,£) be a C°°-function defined on U x U x W 1 , where U is an open subset of IR n . 
Definition 1.1. The function a belongs to the class S™g with p, 5 £ [0, 1] and melif 

(1.1) sup \d^dja{x,y,i)\ < C W7 (1 + iei) m+5(|aWI) - pM 

(i,i/)eJf 

for each compact set C £7 x U and all multi-indices a,/?, 7, where CV ;Qi/ g i7 are some 
positive constants. 

Definition 1.2. An operator A : Cq°(U) h-> V'{U) is said to be a pseudodifferential operator 
of class if 

(c) its Schwartz kernel ^4(x,y) is infinitely differentiable outside the diagonal {a; = £/}, 
(ci) A(x, y) = (2ir)~ n J e l ^ x ~ y '"^a{x, y,£) d£ with some a £ S^j in a neighbourhood of the 
diagonal. 

The function a in (ci) is called an amplitude, and the number m is said to be the order of 
the amplitude a and the corresponding PDO A. Note that for amplitudes of order m > —n 
the integral in (ci) does not converge in the usual sense. However, it is well defined as a 
distribution in x and y. 

Let S~°° := nimGR'S'p^' an d let be the class of operators with infinitely differentiable 
Schwartz kernels. If a £ S~°° (that is, if a and all its derivatives vanish faster than any power 
of |£| as |£| — > 00) then the corresponding PDO A belongs to \l/ _oc . The classical theory of 
PDOs is used to study singularities. Therefore one usually assumes that a is defined modulo 
S~°° and that x is close to y. 

Let a £ S™ s and cij £ for some p,5 £ [0,1], where rrij — > — 00 as j — > 00. We shall 
write 



[1.2) a ~ ' 1^1 ^ 00 



if a — aj £ S™£ where n fc — >■ —00 as — > 00. Such series . aj are called asymptotic. 

If rrij — > —00 then for every collection of amplitudes aj £ S I there exists an amplitude a 
satisfying (11.21) . Obviously, if a' is another amplitude satisfying (II. 2p then a — a 1 £ S~°° (or, 
in other words, (11.21) defines a modulo S~°°). 

Lemma 1.3. Let z T := x + r(y — x) where r £ [0, 1]. If 5 < p and a £ S™ s then 



modulo an infinitely differentiable function, where (ta,t(z,£) is an amplitude of class S™ 5 
given by the asymptotic expansion 

(1.3) ~ E— ffll — W?3fV*.V,0 , leKoo. 

z — ' a! p! y=x=z 

a,/9 
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Sketch of proof. Expand the amplitude a by Taylor's formula at the point (x,y) = (z T ,z T ), 
replace (y — x)e^ x ~ y ^'^ with iV^e 4 ^' - ^ and integrate by parts with respect to £. □ 

The amplitude 071^(2, £) is called the t -symbol of the PDO A. It is uniquely defined by 
the operator A modulo S~°°. The 0-symbol is usually called just the symbol and is denoted 
a a- The |-symbol and 1-symbol are said to be the Weyl and the dual symbol respectively. 

In the theory of PDOs, properties of operators are usually described and results are stated 
in terms of their symbols. The following composition formula plays a key role in the symbolic 
calculus. 

Theorem 1.4. Let A e ty™} and B e #™£. If 5 < p then the composition AB is a PDO of 
class \]>™i +m2 whose symbol admits the asymptotic expansion 

(1.4) OT AB (X,0 ~ ^^^(1,^^(1^), KI-KX). 

a 

Sketch of proof . From the inversion formula for the Fourier transform it follows that the 
Schwartz kernel of AB is given by (ci) with the amplitude a(x,y,£) = aA,(%,£)o-B,i(y,£)- 
Applying Lemma [1.31 with r = to a, we obtain (11. 411 . □ 

Remark 1.5. Theorem 11.41 implies, in particular, that the resolvent of an elliptic differential 
operator is a PDO. Using (II. 4ft . one can also show that a PDO of order m maps Wf H ^0° 
into W2~ m , where are the Sobolev spaces. 

Note that in the above lemmas the condition 5 < p is of crucial importance; if it is not 
fulfilled then the terms in the right hand sides of 111. 3ft and II 1.4ft do not form asymptotic 
series. 

Clearly, the phase function (x — y) ■ £ in (ci) depends on the choice of coordinates on U. 
Passing to new coordinates x and y, we obtain 

A(x,y) = (2tt)-" J e< x ^~y^<a(x(x)My),0^. 

In a sufficiently small neighbourhood of the diagonal {x = y}, the new phase function 
<f(x, y, ^) = (x(x) — y(y)) ■ £ can be written in the form 

f(x,y,0 = (x-y)- $(x,y)£, 

where $>(x,y) is a smooth n x n-matrix function such that det<&(x,y) 7^ 0. Changing 
variables r\ = $(£, y) we see that 

A(x,y) = (2tt)-" J e^-y^a(x,y, V )d V , 

where 

a(x,y,v) = |det$(x,y)|~ 1 a(x(x),y(y),$ _1 (x,y)?7) 

is a new amplitude. Thus Definition 11.21 does not depend on the choice of coordinates. 
However, there are two obvious problems. 
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Problem 1.6. If a G S™ 5 then, generally speaking, the new amplitude a belongs only to the 
class S™ S i with 5' := max{<5, 1 — p}. If p < | then 5' > p and the above lemmas fail. Thus 
for 5 < p < \ it is impossible to define PDOs of class S™ s on a manifold and to develop a 
symbolic calculus using local coordinates. 

Problem 1.7. If max{5, 1 — p] < p then the "main part" of the symbol a a (called the 
'principal symbol of A) behaves as a function on the cotangent bundle under change of coor- 
dinates. However, lower order terms in (jl.3p do not have an invariant meaning. Therefore, 
the coordinate approach does not allow one to study the subtle properties of PDOs, which 
depend on the lower order terms. 



2. Linear connections 

The above problems do not arise if we define the phase function (x — y) ■ £ in an invariant 
way, without using local coordinates. It is possible, in particular, when the manifold is 
equipped with a linear connection. In this section we shall briefly recall some relevant 
definitions and results from differential geometry. 

Let M be an n-dimensional C°°-manifold. Further on we shall denote the points of M 
by x, y or z. The same letters will be used for local coordinates on M. Similarly, £, rj and 
C will denote points of (or the dual coordinates) on the fibres T*M, T*M and T*M of the 
cotangent bundle T*M. 

We are going to consider operators acting in the spaces of K-densities on M, /t£l Recall 
that a complex-valued "function" u on M is said to be a K-density if it behaves under change 
of coordinates in the following way 

u(y) = | det{dx l /dy j }\ K u(x(y)) . 

The usual functions on M are O-densities. The ^-densities are sections of some complex 
linear bundle f2 K over M. We denote by C°°(M; f2 K ) and C^°(M; fi K ) the spaces of smooth 
^-densities and smooth K-densities with compact supports respectively. If u G Cq°(M; f2 K ) 
and v G C°°(M;f2 1_K ) then the product uv is a density and the integral J M uvdx is 
independent of the choice of coordinates. This allows one to define the inner product (u, v) = 
J M uvdx on the space of half-densities C^°(M; Vt 1 / 2 ) and to introduce the Hilbert space 
L 2 (M; fi 1/2 ) in the standard way. 

In this and the next sections we shall be assuming that the manifold M is provided with a 
linear connection V (which may be non-complete). This means that, for each local coordinate 
system, we have fixed a set of smooth "functions" Y l - k {x), i,j,k = l,...,n, which behave 
under change of coordinates in the following way, 

l p,q 

The "functions" T l - k {x) are called the Christoffel symbols. They can be chosen in an arbitrary 
way (provided that (12. ip holds), and every set of Christoffel symbols determines a linear 
connection of M. 
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A linear connection T is uniquely characterized by the torsion tensor T % - k := Y % - k — n . and 
the curvature tensor 

Rjki := ®v k T]j — dy'^lj + r^r^ ~~ r^- . 

p p 

If both these tensors vanish on an open set U C M then one can choose local coordinates 
on a neighbourhood of each point x G U in such a way that Y l - k = 0. Such connections are 
called flat. A connection T is called symmetric if T l - k = 0. 

Let z/ = v k (y) d y k be a vector field on M . The equality (12. ip implies that 

(2-2) V, := XX(y)^ + >2 I /,!//) !//)'/; ^, 

ft 

is a correctly defined vector field on T*M. The vector field (12. 2 p is called the horizontal lift 
of ^ . The horizontal lifts generate a dimensional subbundle HT*M of the tangent bundle 
TT*M over T*M, which is called the horizontal distribution. The vertical vector fields 
d m , . . . , d Vn generate another n-dimensional subbundle VT*M C TT*M which is called 
the vertical distribution. Since HT*M n VT*M = {0} , the tangent space Tr y i j\T*M at 
each point (?/, 77) G T*M coincides with the sum of its horizontal and vertical subspaces. 
Obviously, the horizontal subspaces depend on the choice of V whereas the vertical subspaces 
do not. 

A curve in the cotangent bundle T*M is said to be horizontal (or vertical) if its tangent 
vectors belong to HT*M (or VT*M). For any given curve y(t) C M and covector rjo G 
T*,q,M there exists a unique horizontal curve (y(t), n(t)) C T*M starting at the point 
(y(0), r)o) . It is defined in local coordinates y by the equations 

j fVj (t) - Y, T Uy(t))y\t)m{t) = o, Vj = i,...,n, 

i,k 

and is called the horizontal lift of y(t) . The corresponding linear transformation 7] — >■ 77 (t) 
is said to be the parallel displacement along the curve y(t). By duality, horizontal curves 
and parallel displacements are defined in the tangent bundle TM (and then in all the tensor 
bundles over M). 

A curve y(t) C M is said to be a geodesic if the curve (y(t),y(t)) C TM is horizontal or, 
equivalently, if 

y\t) + X)rj(y(t))y'(t)^(t) = 0, V* = l,...,n. 

in any local coordinate system. If C/ x is a sufficiently small neighbourhood of x then for 
every y <E U x there exists a unique geodesic jy tX (t) such that 7^(0) = x and 7^(1) = y. 
The mapping U x 3 y 1— > ^^(O) G T X M is a bijection between ?7 X and a neighbourhood of the 
origin in T X M, and the corresponding coordinates on U x are called the normal coordinates. 
In the normal coordinates y centred at x we have 7 3/ ,z(£) = x + t(y — x), so that A fy, x {t) — y — x 
for all t G [0, 1]. 

Let : T a *M — >■ T*M be the parallel displacement along the geodesic y y>x , and let 
Yy,x = Idet^ajj. One can easily check that T y x is a density in y and a (— l)-density in 
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x (the map w x — > T VjX w x is the parallel displacement along j VjX between the fibres of the 
bundle Q). Note that $> yiX and T y>x depend on the torsion tensor, whereas the geodesies are 
determined only by the symmetric part of T. 

Given local coordinates x = {x 1 , . . . , x n }, let us denote by V, the horizontal lifts of vector 
fields d x i. For a multi-index a with |a| = q, let V" = ^ ^ . . . Vi 8 where the sum is 
taken over all ordered collections of indices i\, . . . ,i q corresponding to the multi-index a. 
The following simple lemma can be found, for instance, in |Sa2} Section 3]. 

Lemma 2.1. If a G C°°(T*M) then a(y,^ VtX ^) admits the following asymptotic expansion, 
(2.3) a(y,$„.*0 ~ J2^y,x V >( x ^)i V^*- 

a 

Sketch of proof. Write down the left hand side in normal coordinates y centred at x and 
apply Taylor's formula. □ 

Remark 2.2. If a is a function on T*M then the "hypermatrix" {d^V x a(x, £)}\ a \ = q ; \p\ =p 
behaves as a (p, g)-tensor under change of coordinates. Therefore all the formulae in the 
next section have an invariant meaning and do not depend on the choice of coordinates. 

3. PDOS: A COORDINATE-FREE APPROACH 

Definition 3.1. We shall say that an amplitude a defined on M x T*M belongs to the class 
S™ S (T) with p, 5 G [0, 1] and m G R if 

(3.1) sup \d a x dN h V l2 ...V lq a{x,z,Q\ ^ C K ., a>p , q (1 + |£|) m+ *H+«>-'M 
(x,z)eK 

for each compact set K C M x M, all multi-indices a, (3 and all sets of indices i\,...,i q , 
where are horizontal lifts of the vector fields d Zk and CK- a ,/3,-y are some positive constants. 

From the definition of the horizontal lifts it follows that a G S™ S (T) with 5 ^ 1 — p if and 
only if a satisfies fll.ip in any local coordinate system. In this case the class S™ S (T) is the 
same for all linear connections T. If 8 < 1 — p then S™ S (T) depends on the choice of T. Note 
that fll.lft is a particular case of (13. ip . in which the connection T is flat. 

Let us fix a sufficiently small neighbourhood V of the diagonal inMxM and define 
z T = z T (x,y) = jy, x (T) , where r G [0, 1] is regarded as a parameter. Consider the phase 
function 

(3.2) <p T ( x , C, y) = -<%,.(t),C}, (x,y)eV, c, e v: .\/. 

where (•, •) denotes the standard pairing between vectors and covectors. The function ip T 
is invariantly defined and, by the above, coincides with (x — y) ■ ( in normal coordinates y 
centred at x. 

Definition 3.2. An operator A acting in the space of ^-densities on M is said to be a PDO 
of class ^(fi K ,r) if 

(c) its Schwartz kernel A(x, y) is infinitely differentiable outside the diagonal {x = y}, 
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(c 2 ) A(x, y) = (2tt) n p K , T (x,y) J T „ M e Vf>T ^ ,v ' a{z s \ z T , () d( in a neighbourhood of the 

diagonal, where a G S™ S (T), p K>r := Tj^T"^ and s,t G [0,1] are some fixed 
numbers. 

Remark 3.3. If y are normal coordinates centred at x then <p T (x,(,y) = — y) • C an d 
the integral (C2) takes the form (ci). However, Definition 11.21 assumes that x and y are the 
same local coordinates on U, whereas the above identity holds if we choose coordinates y 
depending on the point x. 

Remark 3.4. The weight factor p K)T is introduced for the following two reasons. 

(1) It makes the definition independent of the choice of coordinates ( in the cotangent 
space T* t M. 

(2) Because of this factor, the Schwartz kernel behaves as a (1 — ft)-density in y and 
ft-density in x, that is, (c 2 ) defines an operator in the space of K-densities for all 
k G K and all s, r G [0,1]. In particular, this allows us to consider PDOs in the 
Hilbert space L 2 (M, Vt 1 ^ 2 ) and to introduce Weyl symbols (corresponding to r = ~). 

One can replace p KiT in Definition 13.21 with any other smooth weight factor p(x, y) which 
behaves in a similar way under change of coordinates. The precise choice of the weight factor 
seems to be of little importance, since all formulae in the symbolic calculi corresponding to 
different weight factors p and p can easily be deduced from each other by expanding the 
function p~ 1 p into an asymptotic series of the form (I2.3p . replacing jy >x (z T ) e 1 ^^^'^ with 
iV^ e vpT ^ x ' < '' v ' and integrating by parts with respect to £. 

Lemma 3.5. If 5 < p and a G S™ S (T) then for all s, r G [0, 1] 



and 



Pk,t I e^<,y) a {z s ;z T ,()d( = p K , T / e^<,y)<j AT ( ZT X)d( 

JT* M JT* M 

Zf Zt 

p K , T [ e^^a4, T (z T ,C)dC = p K , s [ e*f"M>ri<T Ata (z s1 Qd{ 

JT* M JT* M 

Z T Z S 

modulo C°° -densities, where o~a,t cmd cta,s are amplitudes of class S™ S (T) given by the as- 
ymptotic expansions 

a 

a A , s (x,0 ~ ^ (-') W 0j -«)'"' ffv^fofl, |£| ->oo. 

a 

Sketch of proof . The first identity is proved by applying ( 12.31) with x = z T and y = z s to the 
function a(-; z T , £) with fixed (z T , (), substituting 7« S)Zt e ltfT — (r — s) V^e tLpT and integrating 
by parts. The second is obtained in a similar way, after changing variables ( = Q Zs ,z T C ■ D 

Lemma [3.51 shows that Definition 13.21 does not depend on the choice of r and s, and that 
every PDO A is defined modulo vl/ -00 by its r-symbol o~a,t- The other way round, for each 
linear connection T, the r-symbol o~a,t is determined by the operator A modulo S~°°. 
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If A G ^™ S (Q K , T) then, in a similar way, one can show that 

(3.3) * A .A*,Z) ~ E i — d^a AiT (x,0, |e|->oo, 

z — ' a! 

where A* is the adjoint operator acting in the space of (1 — ft)-densities. In particular, for 
the Weyl symbols we have 0a*,i/2 — va,i/2 £ S* 00 for all kg8. 

Remark 3.6. The full r-symbol o"a, t depends on T and r. If max{<5, 1 — p] < p then all the 
r-symbols o~a,t corresponding to different connections T coincide with the principal symbol 
of A modulo a lower order term. However, in the general case it seems to be impossible to 
define a principal symbol of A without introducing an additional structure on the manifold 
Mora global phase function (see Subsection I6.8j) . 

Let T K (x, y, z) := Tj"" T 2 ~ X K T£-« , ^(x, £; y, z) := - - (%,,, and 

\p\<\p\ v 

where y and z are normal coordinates centred at x. The functions P^l G C°°(T*M) are 

polynomials in £ ; we shall denote their degrees by dp . 

One can easily show that = 0, Pi^j = and di K ^ < min{ |7| } for any connection 

T. Moreover, if V is symmetric then dp < min{ I7I, (|/3| + |7|)/3 } |Sa2| Lemma 8.1]. 

Theorem 3.7. Let A G ^™}(n K ,T) and B G ^f(fi K ,T) , w/jere p > <5. Assume, 2n 
addition, that 

(1) either p> 1/2, 

(2) or i/ie connection T is symmetric and p > 1/3 , 

(3) or at /east one of the PDOs A and B belongs to ^ (Q K ,T). 

Then AB G ^} +m2 {Q K ,T) and 

(3.4) a AB (x,0 ~ E A ^(*.0 £>f %M ^V^(x,0 , |£| 00. 

a,/3,7 P ' T 

Proof of Theorem 13.71 is similar to that of Theorem 11.41 (see [Sa2, Section 8]). In particular, 
if the connection T is flat then P^l = as + I7I > 1 and (13. 4p turns into (jl.4p . 

Remark 3.8. The conditions on p and the estimates for d^ imply that the terms in the 
right hand side of (13.41) form an asymptotic series. It is plausible that the composition 
formula (13. 4p holds whenever the orders of the terms in the right hand side tend to —00 as 
l a l + \P\ + \l\ ~~ ^ 00 • However, it is not clear how this can be proved. 

Remark 3.9. Coefficients of the polynomials P^l are components of some tensors, which are 
polynomials in the curvature and torsion tensors and their symmetric covariant differentials. 

In the same way as in the local theory of PDOs, Theorem 13.71 implies standard results on 
the boundedness of PDOs in the Sobolev spaces and allows one to construct the resolvent 
of an elliptic operator in the form of a PDO. 
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4. Functions of the Laplacian 

In this section we assume that M is a compact Riemannian manifold without boundary 

and denote |£| x := 9^( x ) where £ G 7 .\/ and {<7 tJ '} is the metric tensor. It is 

well known that there exist a unique symmetric connection T g on M, called the Levi-Civita 
connection, such that the function \£\ x is constant along every horizontal curve in T*M. 
Denote by A the Laplace operator acting in the space of half-densities; in local coordinates 

Au(x) = g K - l (x)^2d xi (g(x)g^(x)d xj (g- K (x)u(x))), 

where g :— | det the canonical Riemannian density. Let v be a self-adjoint first 

order PDO such that —A + v > 0, and let A v := a/ - A + v. The operator A u is a PDO of 
class ty\ o whose symbol coincides with \^\ x modulo S° in any local coordinate system. Thus 
we have A v G (fi 1 / 2 , T) for any linear connection T. 

Definition 4.1. If p G (0, 1], let S™ be the class of infinitely differentiate functions u on 
M. such that 

(4.1) < C,-(l + | S |r-^, Vj=0,l,..., 

where Cfc are some constants. 

A natural conjecture is that the operator u{A v ) is a PDO whenever u G S™. If it is 
true then the symbol of to {Ay) should coincide with w(|£|a.) modulo lower order terms. If 
p < 1/2 then, generally speaking, this function does not belong to the class S™ s with p > 5 
in any local coordinate system. However, since its horizontal derivatives corresponding to 
the Levi-Civita connection are equal to zero, we have w(|£| x ) G S™ S {T S ) for all u G S™. In 
particular, this implies the following 

Lemma 4.2. Let r G [0,1) and U T {t) := exp(itA T u ). Then U T {t) G ^ Ti0 (fi 1/2 , T g ) for all 
t G R and G- UT(t) {x,0 = e u ^* 6 (r) (t,x,£), where o (t) G C°°(R x T*M) and d?b^ G S? >0 /or 
a// = 0, 1, . . . and each fixed t. 

Sketch of proof. Write down U T {t) formally as an integral (C2) with an unknown symbol of 
the form e**'^* b^ T \t, x, £), substitute the integral into the equation d t U T {t) = iAlU T {t), apply 
the composition formula (13 ,4p to Aj,U T {t) and equate terms of the same order in the right 
and left hand sides. □ 

Using Lemma [4. 2\ one can construct other functions of the operator A v . 

Theorem 4.3. If u G S™ then cj{A u ) G ^ {Q^ 2 ,T S ) and 

00 

(4-2) <r u{Au) ~ u(\Z\*) + ^c jtU (x,£)u>V\\£\ x ), Kl^oo, 
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where u^) := d J s uj and Cj„(x,£) G S® . The functions Cj jU are determined recursively by the 
equations 

(4-3) a A u{x,i) = l^ + J^jj^ie^CjA^O- 

i=i ^ J) ' 

Sketch of proof . Define u T (s) = u(s 1 ^ T ), and let u^-{t) be the Fourier transform of tu T . Then 

u(A„) = (2n)- 1 J u£(i) e itAl dt. 

Let G Cq°(R) be equal to 1 in a neighbourhood of the origin and have support contained 
in a small neighbourhood of the origin. Consider the operators 

= (27T)- 1 / q(t) u£(t) e itAl dt, 



w 2 (A,) = (27T)- 1 y (1 - ? (t)) 5£(t) dt. 
By integration by parts, the operator ^(A,) can be written as 

uj 2 {A u ) = (27T)- 1 A; fc / D k t ((1 - ? (t)) ££(*)) dt. 



Since may be chosen arbitrarily large, this shows that ^(A,) has an infinitely smooth 
kernel. By Lemma I4.2[ the operator ui(A v ) is a PDO whose symbol coincides with 



oc 



(27T)- 1 / ^w^e^fc^t^Odt. 



Expanding q(t)b^(t,x,^) by Taylor's formula at t — 0, we see that the symbol of u;i(A„) 
admits an asymptotic expansion of the form f!4.2[) with some functions Cj >u . These functions 
do not depend on u and can be found by substituting u>(s) = s k with k = 1,2,... This leads 
to □ 

Definition 4.4. If p G (0, 1], let ST(g) be the class of C°°-functions on T*M which admit 
asymptotic expansions of the form 

oo 

(4.4) a{x,0 ~ ^ Ci (x,0^i(|^y , Kl^oo, 

i=o 

where Cj G 5° 0) G 5^ with mo = m and — > — oo. Denote by ^/"(fi 1 / 2 , g) the class of 
PDOs acting in the space of half-densities whose r g -symbols belong to >S™(g). 

Theorem 14.31 immediately implies that oj{A v ) G ^^(fi 1 / 2 , g) whenever co G S™. The other 

way round, any PDOs of class ^^(O 1 / 2 , g) can be represented in terms of functions of the 
operator A v . 
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Lemma 4.5. For each A G \l/™(fi 1//2 , g) there exist PDOs C^ v G and functions Uj G Sp 
such that Iq — m, lj — y — oo and 

oo 

(4.5) A ~ J2 C ^"M»), 

j=0 

where the ~ sign means that the Schwartz kernel of the difference A — Yl^o^j^^ji^) 
becomes smoother and smoother as k — > oo. 

Sketch of proof . Assume that (14. 4p holds and denote by Co the PDO with symbol Cq(x,^). 
Theorems O and Ed imply that A = C tu (A u ) + A^ where Afi G ^(fi 1/2 ,g) with 
l\ ^ maxlm^mo — p}. The same arguments show that A„ = Ci^Cji(A v ) + A^ where 
c i,v G *i,o> ^i G S p and G ^p( fil/2 >g) where l 2 ^ max{m 2 J 1 - p}. Repeatedly 
applying this procedure, we obtain a sequence of operators Au G ^> l p k (Vl l l 2 , g) such that 

A — Av = Yl'jZo Cj, u u)j(A u ), where Cj )V and ujj satisfy the required conditions and l k — > — oo 
as /c — >■ oo. □ 

Since ^(A^) u; 2 (A,) = Wiu^ (A,) for any two functions U\ and u 2 , combining Theorem 
13.71 and Lemma 14.51 we obtain 

Corollary 4.6. If A G ^^(fi 1 / 2 , g) and 5 G ^>™ 2 (fi 1/2 ,g) t/ien i/ie composition AB is a 
PDO of class \l/£ ll+m2 (f2 1 / 2 ,g) whose symbol admits the asymptotic expansion (13 ,4p . 

Remark 4.7. Under the conditions of Corollary 14. 6[ the estimates on <ig i, obtained in Section 
[3] do not directly imply that H3.4[) is an asymptotic series, as it seems to contain terms 
of growing orders. However, these "bad" terms cancel out due to the symmetries of the 
curvature tensor. It would be interesting to find a direct proof of Corollary 14.61 which does 
not use Lemma 14.51 (a relevant problem was mentioned in Remark 13. 8p . 

From the above results it follows that the restriction of the operator u(A v ) to an open 
subset of M is determined modulo \]> -00 by the restrictions of the metric g and the operator 
v to this subset. More precisely, we have the following 

Corollary 4.8. Let v G C™(M), and let {v} be the corresponding multiplication operator. 
Consider the operator A„ generated by another metric g and another first order PDO v. If 
g = g on the support of the function v and v {v} = v {v} then 

{v} (u(A v ) - u(A p )) G and (u(A„) - u(A p )) {v} G 

for every u G S™. 

Sketch of proof. The multiplication operator {v} is a PDO with symbol v(x), which belongs 
to ^(fi 1 / 2 , g). Applying Lemma I4T51 and Corollary 14. 6[ we see that {v}(co(A u ) — u(A„)) 
and (u(A v ) — uj(Ai,)){v} are PDOs whose full symbols are identically equal to zero. □ 

Remark 4.9. In a similar way, it is possible to define the classes ty™(Q K ,g) which consist of 
PDOs acting in the space of ^-densities. Theorem 13.71 implies that A G ty™(fl, K , g) if and 
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only if g 1 / 2 K Ag K x l 2 G \l/™(f2 1//2 , g). Using this observation, one can easily reformulate all 
results of this section for operators A G \I/™(f2 K ,g). 

5. An approximate spectral projection 

In applications, one often has to deal with functions of an operator which depend on 
additional parameters. It is more or less clear that the results of the previous section can be 
extended to parameter-dependent functions u under the assumption that the estimates (14 .1}) 
hold uniformly with respect to the parameters. Therefore, instead of formulating general 
statements, we shall consider an example which is of particular interest for spectral theory. 

Further on we assume that A > and denote by \E' _00 (A) the class of parameter-dependent 
operators with infinitely smooth Schwartz kernels A\(x,y) such that 

lim \ p \d^d^A x (x,y)\ = 

for all multi-indices a, (3 and all p = 1, 2, . . . Similarly, let S~°°(\) be the class of parameter- 
dependent amplitudes a\ such that 

{X + W sup |^V£V> A (y;x,£)| asA+|£|^oo 

for all multi-indices a, (3, 7 and all p — 1, 2, . . . 

Let us fix a small e > and a nonincreasing function / G C°°(R) such that 

1 if s € 0; 



^ ~" 1 if s ^e, 

and ^ f(s) < 1 for all s£i If p G (0, 1] and A > 0, let 

Xp (\,s) := /(A-'(a-A)). 

For each fixed A > 0, the function \ p vanishes on the interval [A+eA p , 00), is identically equal 
to 1 on the interval (—00, A] and smoothly descends from 1 to on the interval [A, A + e\ p \. 
Since this functions differs from the characteristic function of the interval [—00, A] only 
on the relatively small interval (A, A + e\ p ), the operator x(\A v ) can be thought of as 
an approximate spectral projection of A v corresponding to (—00, A]. The standard elliptic 
regularity theorem implies that the operator x(A, A u ) has an infinitely differentiable Schwartz 
kernel for each fixed A. 

The derivatives d{xp{\ s ) are equal to zero outside the interval (A, A + e\ p ). Therefore 

(5.1) \dlx P {\s)\ < ^.(H+A)-", Vj = 0,l,..., 

for all s G M and all A > 1, where Cj are some constants independent of A and s. The same 
arguments as in the proof of Theorem 14.31 show that x(A, A v ) is a parameter-dependent PDO 
whose symbol admits the asymptotic expansion 

00 

(5.2) a x{XAv) ~ x(A,|eU) + ^c^(x,Ox ' } (A,|eU), A + |^| ^ 00, 

3=1 

where Cj v (x,£) are the same functions as in (14.31) and x^ denotes jth s-derivative of the 
function x- 
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Note that the functions X (A, |£U) belong to S~°° for each fixed A. However, their rate 
of decay depends on A. The asymptotic expansion (15. 2p is uniform with respect to A; it 
defines c x a,A v ) modulo S~°°(X). Substituting the terms from ( 15. 2 p into the integral (c 2 ), we 
obtain an asymptotic expansion of the Schwartz kernel of x(A, A u ) into a series of infinitely 
smooth half-densities, which decay more and more rapidly as A — > oo. This expansion defines 
x(X,A„) modulo *-°°(A). 

Straightforward analysis of the proof of Theorem 13.71 shows that it remains valid in the 
case where one of the operators belongs to ^ and the other is a parameter-dependent PDO 
whose symbol admits an asymptotic expansion of the form ( 15. 2p . In this case ( 13. 4p gives an 
expansion of o~ab as A + |£| — > oo and defines the symbol modulo S~°°(X). 

Now, in the same way as in Section HJ one can show that the composition of parameter- 
dependent PDOs whose symbols admit asymptotic expansions of the form 

oo 

(5.3) a{x,£) ~ ^foOx^KW, A+|e|^oo, Cj G < 

3=0 

is also a parameter-dependent PDO whose symbol is given by (13 .4p modulo S~°°(X). 

Let H u (\) be the spectral projection of the operator A v corresponding to the interval 
(— oo, A). The above results imply the following 

Theorem 5.1. Let v G C^°(M), and let {v} be the corresponding multiplication operator. 
Consider the spectral projections n„(A) and IIp(A) generated by different metrics g, g and 
different first order PDOs u, v satisfying the conditions of Section HI // g = g on the 
support of the function v then 

(5.4) n„(A) {v} (I - %(X + cX p )) G *-°°(A), Vc,p>0. 

Sketch of proof Assume that p G (0,1], and let x(A, s) be defined as above with some 
e < c/3. Then x(A, s) = 1 for s ^ A and x(A, s - e\ p ) = for s ^ A + c\p. It follows 
that 

n i/ (A)x(A,A 1/ )=n„(A) and X {\A -eWI) [I - + c\ p )) =0. 
Consequently, we have 

(5.5) n v (A)) {v} (I - U (X + cX p )) 

= n„(A) x(A, A v )) {v} (I - X (A, Ap - e\ p I)) (I - U,(X + cX p )) 

Since x(A, s) = x(A, s) x(A, s — eX p ), the composition formula implies that 

X (A,A i/ ))M(/-x(A,A p - £ A' ) /)) G vi/— (A). 

Now the required result follows from (15. 5p and the fact that the Schwartz kernel of the 
spectral projection is polynomially bounded in A with all its derivatives (see, for instance, 
[5V1 Section 1.8]). □ 

Remark 5.2. It is not surprising that the operator in the left hand side of ( 15 ,4p has a lower 
order than the spectral projections themselves as A — > oo. However, one would expect its 
norm to decay as a fixed negative power of A, since the perturbation A u — A„ is a more or 
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less arbitrary PDO of order zero. We do not know whether ( 15. 4ft can be obtained by other 
techniques (including that of Fourier integral operators). 

Remark 5.3. All results of this section can easily be extended to a noncompact closed mani- 
fold M. In this case all the asymptotic expansions are unform on compact subsets of M and 
M x M. 

6. Other known results and possible developments 

6.1. Other definitions for scalar PDOs. If T is a linear connection, then the corre- 
sponding symbol of a PDO A can easily be recovered from the asymptotic expansion of 
A(e lipT ( x ^' y ^x( x , y)) as £ — )• oo, where ip T is defined by (13. 2p and x is a smooth cut-off func- 
tion or ^-density (we suppose that x is fixed and that the operator acts in the variable y). 
After that, all the standard formulae of the local theory of PDOs can be rewritten in terms 
of their T-symbols. Moreover, making appropriate assumptions about the asymptotic be- 
haviour of A(e lipT<yX ' c - , ' s, x{ x i y))i one can try to define various classes of PDOs associated with 
the linear connection Y. 

This approach was introduced and developed by Harold Widom and Lance Drager (see 
|Wilj . [Wi2] and |Drj ) . Its main disadvantage is the absence of an explicit formula represent- 
ing the Schwartz kernel of a PDO via its symbol. As a consequence, one has to assume that 
PDOs and the corresponding classes of amplitudes are defined in local coordinates, which 
makes it impossible to extend the definition to p < max{<5, 1 — p}. 

In |Pf 1] . Markus Pflaum defined a PDO in the space of functions by the formula 

(6.1) Au(x) = (2tt)-*/ / X (x,y)e iMx 't' y) a(x,0u(y)dyd£, 

Jt*m Jt x m 

where a(x, £) is a function on T*M of class S™ s , y are normal coordinates centred at x 
and x is a smooth cut-off function vanishing outside a neighbourhood of the diagonal. He 
obtained asymptotic expansions for the symbols of the adjoint operator and the composition 
of PDOs and, in the later paper |Pf2] . extended them to r-symbols. However, the results in 
|Pflt IPf2] are stated and proved with the use of local coordinates and, therefore, the author 
had to assume that max{<5, 1 — p} < p. 

Recall that under this condition the standard results of the local theory of PDOs hold, 
and the only advantage of a coordinate-free calculus is that it helps to fight Problem 11.71 A 
typical example, considered in |Pf 1] . is the PDO with a symbol of the form (1 + |£| 2 ) fc ^ where 
b(x) is a smooth function on M. Formally speaking, this PDO belongs only to the class S™ s 
with m = sup x b(x) and any 5 G (0, 1). But its properties are determined by the values of 
the function b at all points x G M; in a sense, this operator has a variable order depending 
on x G M. In such a situation, it is not sufficient to consider only the principal symbol. One 
has to define a full symbol which can be done with the use of a linear connection. 

It is clear that ( 16. ip differs from Definition l3.2l only by the choice of the weight factor p K T . 
Applying the procedure described in Remark 13.41 one can easily show that 



(6.2) 



°A{x,i) ~ ^P a (x)<9 5 Q a(x,0 , |f|->oo, 

a 
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where a(x,£) is the symbol appearing in (16. ip and P a are components of some tensor fields. 
Using (16. 2p . one can rewrite all the results obtained in |Sa2] in terms of symbols defined 
by (16. ip . This shows that Pflaum's formulae can be reformulated in terms of the horizontal 
derivatives V" and thus extended to the classes ^ r ^ 5 {VL K , Y) and r-symbols. 

In particular, Pflaum's composition formula can be written in the form (13.41) with some 
other polynomials Pp K l- For operators acting in the space of functions and r = 0, this result 
was established by Vladimir Sharafutdinov in [SI]. He chose to give a direct proof instead of 
deducing the formula from (13. 4 j) and (16. 2 j) and, for some reason, considered only the classes 
^i^o- Sharafutdinov gave an alternative description of the polynomials PjQ which may be 
useful for obtaining more explicit composition formulae (this investigation was continued in 
|Ga] ) . He also proved an analogue of (I3.3P in the case k — 1/2 and r = [Si] Theorem 6.1]. 

Remark 6.1. From (16.21) it easily follows that the degrees of the polynomials P^l admit the 
same estimates as dp (see Section [3]). 

6.2. Operators on sections of vector bundles. In |FK1 IPf2l IS2l IWi2j the authors con- 
sidered PDOs acting between spaces of sections of vector bundles over M. In this case, in 
order to construct a global symbolic calculus, it is sufficient to define parallel displacement 
and horizontal curves in the induced bundles over T*M. This can be achieved by introducing 
linear connections on M and the vector bundles over M. After that the results are stated 
and proved in the same way as in the scalar case (further details and references can be found 
in the above papers). 

A more radical approach was proposed by Cyril Levy in [Lej . He noticed that in order 
to develop an intrinsic calculus of PDOs one actually needs only an exponential map, which 
does not have to be associated with a linear connection. In his paper Levy assumed that 
the manifold M is noncompact and is provided with a global exponential map (that is, M is 
a manifold with linearization in the sense of [Bo]). He then defined associated maps in the 
induced vectors bundles and constructed a global coordinate-free symbolic calculus. 

Remark 6.2. All the papers mentioned in this subsection dealt only with symbols whose 
restriction to compact subsets of M belong to S™ s with p > max{<5, 1 — p}. It should be 
possible to extend their results to p < 1/2, using the technique outlined in Section 

6.3. Noncompact manifolds. In order to study global properties of PDOs on a noncom- 
pact manifold M, one has to assume that all estimates for symbols and their derivatives 
hold uniformly for all x G M (rather than only on compact subsets of M, as in Definitions 
11.11 and 13. ip . In |Baj , Frank Baldus defined classes of symbols and developed an intrinsic 
calculus of PDOs on a noncompact manifold M under the assumption that M has an atlas 
satisfying certain global conditions. The statements and proofs in |Baj were given in terms 
of local coordinates, and global results were obtained by considering the transition maps 
between coordinates charts. It is quite possible that these results can be simplified or/and 
improved under the assumption that M has a global exponential map (as in [Lej). 

6.4. Other symbol classes. The paper [Ba] dealt with the more general classes of symbols 
S(m,g) instead of S™ s . The classes S(m,g) were introduced by L. Hormander in |Hlj (see 
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also [H2J). They are defined with the use of coordinates, and in each coordinate system S™ s 
is a particular case of S(m,g). It would be interesting to construct similar classes S(m,g) 
associated with a linear connection (or an exponential map) and to study the corresponding 
classes of symbols and PDOs. 

Remark 6.3. Note that the introduction of "coordinate" classes S(m,g) does not help to 
resolve Problem II. 61 The relation between these "coordinate" classes and the classes S™ S (T) 
was discussed in [Sa2| Remark 3.5]. 

6.5. Operators generated by vector fields. Let v := {ux, v-i, ■ ■ ■ , u n } be a family of 
smooth vector fields Uj on M which span T X M at every point x G M. Consider the corre- 
sponding first order differential operators d Uj and denote 

where q = \a\ and the sum is taken over all ordered sets of indices ji, . . . ,j q corresponding to 
the multi-index a = (ax, • • • , «„). In other words, can be thought of as the symmetrized 
composition of d v . . 

The family v generates a unique curvature- free connection T u , with respect to which all 
covariant derivatives of the vector fields Uj are identically equal to zero. The r^-symbol of 
d% coincides with a" 1 . . .cr" n , where = ^(^,0 : = (^k,0 (see |Sa2| Example 5.4]). Since 
the functions are constant along horizontal curves in T*M generated by the connection 
T v , the operators d£ and their linear combinations can be regarded as constant coefficient 
operators relative to the connection T v (or to the family of the vector fields v). 

This observation was used by Eugene Shargorodsky in |Sha] . where he developed a com- 
plete theory of pseudodifferential operators generated by a family of vector fields v. He 
introduced anisotropic analogues of classes S™ s , proved the composition formula for the cor- 
responding classes of PDOs, defined semi-elliptic operators associated with the family u, and 
constructed their resolvents. All the results in |Sha] were obtained for operators acting on 
sections of vector bundles equipped with linear connections (see Section [672]) . 

6.6. Operators on Lie groups. In [RTJ, the authors defined full symbols of scalar PDOs 
on a compact Lie group M in terms of its irreducible representations and developed a calculus 
for such symbols. It would be interesting to compare their formulae with those obtained by 
introducing an invariant linear connection T on M and applying the methods of |Sa2] or 
EM- 

6.7. Geometric aspects and physical applications. The importance of intrinsic ap- 
proach in the theory of PDOs for quantum mechanics is explained in the excellent review 
|Fuj by Stephen Fulling. Further discussions can be found in the PhD thesis |Guj . Various 
geometric applications are considered in [BNPWJ and |Voj . We refer the interested reader 
to the above papers and references therein. 

6.8. Global phase functions. It is worth noticing that one does not need a linear connec- 
tion or even an exponential map to define PDOs on a manifold in a coordinate-free manner. 
It is sufficient to fix a globally defined phase function satisfying certain conditions. 
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Namely, let ip(x;y,T]) be an infinitely differentiable function on M x T*M such that 

Im <p(x; y,rj)^0, tp(x\ y, Xr]) = A ip(x; y, rj) 

for all x G M, (y, 77) G T*M and A > 0, and 

<p(x;y,ij) = (x - y) ■ rj + 0(\x - y\ 2 \v\) , x -> y , 

in any local coordinate system. If a(x; y, if) is a smooth function on M x T*M such that 
a G S™ in any local coordinate system then 

A{x,y) := J e iv{x ' y ' n) a(x;y,ri)dr] 

is the Schwartz kernel of a PDO A G acting in the space of functions. Moreover, there 
exists an amplitude a v (y,rj) independent of x such that 

A(x,y) - J e^y^a^y,^)^ G C°°{M x M) , 

and this amplitude a v is uniquely defined by A modulo S~°°. The operator A belongs to 
V&To if and only if G S™ in any local coordinate system. 

Remark 6.4. For a real-valued phase function f these are standard results of the theory 
of Fourier integral operators (see, for instance, |Shuj Section 19]). Complex- valued phase 
functions were considered in |LS Vj . 

It is natural to call a v the ip-symbol of the operator A. Clearly, all the standard results 
of the classical theory of PDOs can be rewritten in terms of their (^-symbols. In particular, 
if A, B G V&To then the ^-symbol of the composition AB is determined modulo S~°° by an 
asymptotic series which involves ^-symbols of A and B and their derivatives. Similarly, the 
if-symbol of the adjoint operator A* is given by a series involving the derivatives of (^-symbol 
of A. 

Obviously, the same formulae remain valid under milder assumptions about the symbols. 
Thus it should be possible to introduce symbol classes associated with the phase function cp 
and develop a symbolic calculus in these classes (as was done in [Sa2] for the special phase 
function (p T generated by a linear connection). 

Such a general approach may allow one to extend results of Section [5] to other elliptic 
operators. It may also be useful for the study of solutions of hypoelliptic equations and 
operators on noncompact manifolds. 
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